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Nature of the phase transition of regularly frustrated vector spin systems in three dimensions is 
investigated based on a Ginzburg-Landau-type effective Hamiltonian. On the basis of the variational 
analysis of this model, Onoda et al recently suggested the possible occurrence of a chiral phase, 
where the vector chirality exhibits a long-range order without the long-range order of the spin 
[Phys. Rev. Lett. 99, 027206 (2007)]. In the present paper, we elaborate their analysis by 
considering the possibility of a first-order transition which was not taken into account in their 
analysis. We find that the first-order transition indeed occurs within the variational approximation, 
which significantly reduces the stability range of the chiral phase, while the chiral phase still persists 
in a restricted parameter range. Then, we perform an extensive Monte Carlo simulation focusing on 
such a parameter range. Contrary to the variational result, however, we do not find any evidence 
of the chiral phase. The range of the chiral phase, if any, is estimated to be less than 0.1% in the 
temperature width. 
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I. INTRODUCTION 

In vector spin systems, frustrations often induces non- 
collinear or noncoplanar spin structures. Such canted 
spin structures generally accompany the order of the chi- 
rality [lj. Two types of chirality has been discussed 
in the literature. One is a scalar chirality which is 
denned as a scalar product of three Heisenberg spins, 
X ~ Si ■ (Sj x Sk)- The scalar chirality takes a nonzero 
value for a noncoplanar spin configuration. By contrast, 
a vector chirality, which is a target of this paper, is de- 
fined as a vector product of two Heisenberg (or XY) 
spins, k ~ Si x Sj. It takes a nonzero value even for a 
noncollinear but coplanar spin configuration. The order- 
ing of the vector chirality is realized in, e.g., conventional 
helical magnets. 

Although the chiral order inevitably appears in non- 
collinear or noncoplanar spin ordered states, it can be re- 
alized in principle without accompanying the long range 
order of the spin. In such chiral ordered but spin disor- 
dered state, spin correlation lengths are kept finite while 
the chirality shows a long-range order. 

In the past, the existence of such a chiral phase has 
been discussed for several frustrated vector spin systems 
including spin glasses 0-Q and regularly frustrated mag- 
nets [q-UOf . For example, it has been suggested in the 
three-dimensional (3D) Heisenberg spin glass that the 
glass-order of the scalar chirality takes place at a tem- 
perature higher than that of the spin-glass order 
In regular systems, it has been suggested that the two- 
dimensional (2D) fully frustrated XY models exhibits the 
ordering of the vector chirality at a temperature higher 
than that of the spin Kosterlitz-Thouless transition p-Q ■ 
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Such occurrence of separate chiral and spin transitions is 
often called "spin-chirality decoupling" . 

In regularly frustrated 3D systems, however, there has 
been no clear evidence of such an intermediate chiral 
phase so far. For classical Heisenberg and XY antifcr- 
romagnets on the 3D stacked-triangular lattice, Monte 
Carlo (MC) simulations suggested the occurrence of a 
single magnetic phase transition from a paramagnetic 
phase to a helical magnetic phase jllT - fla ]. Based on 
a renormalization-group (RG) analysis, Kawamura sug- 
gested that the phase transition of noncollinear magnets 
could belong to a new "chiral" universality class distinct 
from the standard O(N) Wilson-Fisher universality class, 
whereas the transition could also be of first-order depend- 
ing on the parameter values of the system pjj. Some 
supports to this scenario were reported from field theo- 
retical approaches [H,[l9| and Monte Carlo (MC) simula- 
tions |L2h141 Il9| . Experimental measurements of relevant 
critical exponents also seem consistent with such a chi- 
ral universality class [20-24]. By contrast, some authors 
argued that the transition might be weakly first-order 
[iH liH . |2o| . In any case, though the nature of the tran- 
sition has been somewhat controversial, it has been be- 
lieved that an intermediate chiral phase does not appear 
in 3D regular systems. 

Recently, Onoda and Nagaosa studied the possibility of 
the vector chiral phase in regularly frustrated 3D Heisen- 
berg systems flOl] . Based on a Ginzburg-Landau (GL) 
Hamiltonian describing helical Heisenberg magnets and 
performing variational calculations, these authors sug- 
gested that the Dzyaloshinskii-Moria interaction and/or 
the Coulombic four-spin ring-exchange interaction could 
stabilize the chiral phase even in 3D. 

In this paper, motivated by the recent work by Onoda 
et al, we wish to examine the nature of the phase transi- 
tion of the same GL Hamiltonian as studied by Onoda et 
al by means of a further analytical calculation and a MC 
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simulation on a discretized version of the model, with 
particular attention to the issue of the existence/non- 
existence of an intermediate chiral phase. Note that On- 
oda et al implicit assumed in their analysis a continuous 
nature of the transition, ignoring the possibility of a first- 
order transition jloj]. We see in the present paper that a 
first-order transition indeed occurs within the variational 
approximation, which significantly reduces the stability 
range of the chiral phase. Yet, the variational calculation 
predicts that the chiral phase still persists for a certain 
restricted parameter range. With reference to the results 
of such variational calculation, we also perform extensive 
MC simulations on the lattice discretized version of the 
chiral GL model. In contrast to the variational results, 
MC gives no evidence of the chiral phase. If it exists, the 
stability range of the chiral phase is extremely narrow, its 
width being less than 0.1% in the relative temperature. 

The rest of the paper is organized as follows. In Sec. HH 
we describe the GL model relevant to our present study, 
and briefly review the previous results on the model. In 
Sec. IIII1 we present the results of our variational cal- 
culation taking account of the possibility of a first-order 
transition. In Sec. IIV1 we explain the details of our MC 
simulations. MC results are presented in Sec. |V] for the 
Heisenberg case, and in Sec. I VII for the XY case. Fi- 
nally in Sec. IVIIl we summarize our main results and 
further discuss the possibility of the chiral phase in regu- 
larly frustrated 3D spin systems. Appendices are devoted 
to the details of the variational calculations. 



II. THE CHIRAL GL MODEL 

In this paper, we discus the possibility of the chi- 
ral phase in regularly frustrated 3D vector spin systems 
based on the following Gintzburg-Landau Hamiltonian 

mm, 




+ u(a 2 + b 2 ^ +v (a-tj - a 2 b 



(1) 



where a(r) and b(r) are n-component vector fields associ- 
ated with the noncollinear spin structure at wavevectors 
±Q via 



S(r) = a(r) cos(Q • r) + b{r) sin(Q • r). 



(2) 



In order eq. ([2]) to actually represent the noncollinear spin 
structure, the quartic coupling v should be positive so 
that a and b prefer to be orthogonal to each other. Note 
that in order to bound the free energy, we need to limit 
the range of u and v as 
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FIG. 1: The mean- field phase diagram of the chiral GL model 
in the (u, v) plane. The dotted curve is an expected boundary 
of a first-order transition when fluctuations are introduced. 



This effective Hamiltonian can be derived from a micro- 
scopic spin Hamiltonian with isotropic bilinear interac- 
tions via the Hubbard-Stratonovich transformation [l7j . 
In such a case, the ratio v/u becomes 4/3, while addi- 
tional higher-order terms, which do not explicitly appear 
in (1), are also generated. Hereafter, we call the GL 
Hamiltonian (1) the "chiral GL model". 

In the mean field approximation, a continuous transi- 
tion takes place at r — [17j . When v > 0, the ordered 
phase is a helical magnetic state characterized by 



\a\ 2 



Au 



_L b 



(0 < v < Au). (4) 



When v < 0, the ordered phase is a linearly polarized 
sinusoidal state characterized by 



2u 



(v < 0). 



(5) 



u > and v/u < 4. 



(3) 



When u < or v/u > 4, the free energy is unstable 
and a higher-order term is needed to stabilize it. In such 
a case, the transition generally becomes of first order. 
The mean-field phase diagram is summarized in FigfT] 
It may be worthwhile noting that, although the mean- 
field approximation predicts a continuous transition for 
v/u < 4, fluctuations might change this result leading 
to a first-order transition, especially near the boundary 
v = Au. 

By means of a RG analysis of the chiral GL model, 
Kawamura found a new fixed point distinct from the 
standard Wilson-Fisher O(n) fixed point for certain 
range of the parameters [ll|, I n his analysis, how- 
ever, the possibility of the chiral phase was not consid- 
ered. Since the RG expansion employed in ref.flli Il7j 
was an expansion from dimension four or from the many- 
component limit ri — > oo where the chiral phase is never 
expected to occur, the chiral phase might be missed due 
to an intrinsic limitation of the method employed even if 
it actually exists in 3D in a certain parameter range. 

Other field theoretical approaches supported the exis- 
tence of a new fixed point [l8l . [l9| . They also performed 



3 



a direct Monte Carlo simulation of the chiral GL model 
in case of n — 2 19], and they found evidence of the new 
universality class, although it concerned with a parame- 
ter range different from the target of this paper, which 
corresponded to smaller v/u values. 

Recently, on the basis of a variational approximation, 
Onoda et al predicted that the chiral GL model in 3D 
might exhibit a chiral phase characterized by (a x b) ^ 
with (a) = (b) = 0, if the quartic couplings v and u sat- 
isfy the relation v/u > 4/3 Although there has not 
been clear evidence of the chiral phase in regularly frus- 
trated 3D systems so far, the suggestion by Onoda et al 
promotes us to further examine the possible appearance 
of the chiral phase in the 3D chiral GL model. 



III. ANALYTICAL CONSIDERATION 

In this section, we study the ordering of the 3D chiral 
GL model analytically, either by the variational calcula- 
tion (subsection A) or by the mapping to the nonlinear 
(7 model (subsection B). 



A. Variational Approximation 

In this subsection, following the analysis by Onoda et 
al [id , we study the ordering of the 3D chiral GL model 
based on the variational approximation. Although some 
results were already reported by these authors, we will 
also present them for the sake of completeness. The main 
difference of our analysis from that of Onoda et al is 
that we consider the possibility of a first-order transi- 
tion which was not considered by Onoda et al. In fact, 
a first-order transition is realized within the variational 
approximation, significantly reducing the stability range 
of the chiral phase. 

First, we deal with the case of the Heisenberg spin 
(n = 3), assuming that the macroscopic vector chirality 
appears in the z-direction. The variational Hamiltonian 
for the chiral GL model may be given by 

Ua = lJ dr {( VSa ) 2 + ( VS ^) 2 

+ ry [(Sa x ) 2 + (5a y f + (8b x ) 2 + (8b v ) 2 ] 
+ r ± [(Sa z ) 2 + (Sb z ) 2 ] 

- h K (8a x Sb y - 8a y 8b x )\. (6) 

where da and 8b represent the deviations of the fields 
from their average values with respect to the variational 
Hamiltonian Ho, i.e., 8a(r) = a(r) — (a)o and 8b{r) = 
b(r) — (b)o, where (• • • )o being the average with respect 
to Hq. 

In terms of the new two-component vectors a, /3, 7 



defined by 

a i = ^i a x + b y ), a 2 = -^=(a y - b x ) 

Pl = ^(ay+bx), fo = -j=(-a x + b y ) 

li=a z , l2=b z , (7) 

the variational Hamiltonian Ho can be diagonalized as 
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where V represents the volume of the system. 

Let us denote A = (a)o, B = (/3)o, C = (7)0, with 
8a = a — A, 8/3 — j3 — B, and $7 = 7 — C. Then, the spin 
order parameter and the vector chirality order parameter 
are given by 
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^ V V2C 2 



= A 2 -B 2 



V ^(q 2 +r ll r-hl/A- 



(9) 



B = C = and 



The chiral phase is characterized by A - 

In order to determine the optimal values of the vari- 
ational parameters rii , r±, h K , A, B, C within the present 
variational approximation, we employ the so-called Feyn- 
man inequality, 

T a + (n-n ) >t, (10) 

where T (J-o) is the free energy associated with the 
Hamiltonian H (Ho). The optimal values of the parame- 
ters are then determined by minimizing the l.h.s of (jlOp . 
The detailed form of the l.h.s. is given in Appendix [Al 

Let us assume A\ — m > and other mean values 
are all zero in the ordered state, the vector chirality (a x 
b) pointing along the z-direction. In order to avoid the 
ultraviolet divergence, we introduce here an upper cutoff 
of wavevector A. Various parameters are then rescaled 
as 

T U V 

r= A^' W "2^A' 
_ r_L 
A 2 ' r± - ~ A 2 ' 
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By taking the derivatives of the l.h.s. of (TlCfl) with respect 
to f || , fj_, h K and to, and setting them to zero, we get the 
following conditions for the optimal parameter values 



fx=r + 2(u-~) (m 2 + a 2 a + a 



^ila 2 , 



/^ = 2 ( -5 - 5 ) (to 2 + a 2 - a 2 ) 
= 



2 j // - | ,//- 



//'. 



(13) 
(14) 

(15) 



where the variances of a ~ 7, er 2 ~ cr 2 , are given by 

dq, 



lo q 2 + f || - h„/2 
a 2 =^{6j3 q .6j3- q ) a 



= 2 



/o g 2 + f || + /i K /2 

CT 7 = ' <5 7-q>0 



= 2 



(16) 



/o <7 2 + 5± 

Note that Eq. p3|) has the following two types of solutions 
m = 0, (17) 

_ 2 _ H| _ ^«/ 2 



2(5-5/4)' 



(18) 



The latter case corresponds to the standard helical phase, 
while the chiral phase corresponds to the former case. 

As discussed by Onoda et al, a solution with to = and 
h K ^ is possible if 5 and 5 satisfy a relation 5/5 > 4/3 
[lOj . Setting m = in ([Tip, we get 



h K = h K (35 - 45) / dq 



((? 2 +f||) 2 -^/4 



(19) 



A solution with h K 7^ exists only if 5 is lager than |5. 
The transition to the chiral phase occurs when fu is equal 

(c) 

to f j, satisfying the relation, 



(3v-4u) / dq 



~{q 2 +& ] ) 2 



= 1. 



(20) 



From Eq.(18), one sees that a continuous transition to 
the helical phase occurs at f|| = h K /2. By substituting 



h K = 2f\\ into JT9J), the value of r\ 
to-helical transition is obtained as 

(35-45) / dq- 



at the chiral- 



1 



2fj s) 



= 1. 



(21) 



(12) From Eqs. (l2"0|) and (|2T|). one can show that is al 



ways greater than f \^ . The region f ^ < r || < then 
corresponds to the chiral phase. 

In order to discuss how the ordering proceeds when the 
"temperature" r is varied, we need to calculate the "tran- 
sition temperatures" f c and r s corresponding to f| c ^ and 

( s) (c) ( s) 

fn . Although the inequality f Ji < ffi always holds, 
the relation between the corresponding r c and f s is not 
trivial depending on the values u and v. For a fixed ratio 
5/5, we can show that there is a critical value 5 C such 
that f s > f c for 5 < 5 C and f s < r c for 5 > 5 C (details 
are given in Appendix [A}. Note that only the latter situ- 
ation means the existence of the chiral phase. Indeed, in 
the case of 5 < 5 C , a first-order transition from the para- 
magnetic phase to helical phase occurs directly without 
an intermediate chiral phase. 

Fig. [2] exhibits f c and f s as a function of 5/5 for the 
case of 5 = 0.1 and 0.5. For 5 = 0.1, an inequality 
r c < f s is satisfied for 4/3 < 5/5 < 4 (there may be a 
tiny region of r c > r s in a close vicinity of 5/5 = 4/3 ). 
The chiral phase does not exist for this small value of 5. 
By contrast, for 5 = 0.5, an inequality r c > r s is satisfied 
for 4/3 < 5/5 < 3.68, and the chiral phase exists for this 
large 5. 

So far, we have considered only the case of to = 0. 
However, if a first-order transition really occurs, a con- 
tinuous transition to the chiral phase at f c might be in- 
terrupted by such a first-order transition, and we need 
to examine the case of to 7^ simultaneously, choosing 
the state giving the lower free energy within the varia- 
tional approximation. This point has not been examined 
in Ref. 10]. The first-order transition point is located by 
comparing the free energies (the l.h.s. of (|10p ) of the 
to = and rh ^ solutions. The explicit form of the free 
energy is given in appendix |AJ 

In Figj3l we show the numerically calculated transi- 
tion temperatures as functions of 5 for several values of 
5/5. The first-order transition temperature, ?i, is usually 
greater than r s . It means that the transition to the heli- 
cal phase is of first-order within this variational approx- 
imation. Recall here that, even in the typical 4> A model 
describing a ferromagnet, the same variational approx- 
imation predicts an artificial first-order transition con- 
trary to the reality. Therefore, the first-order nature of 
the transition might also be an artifact of the variational 
approximation employed here. In case of 5/5 = 2.0, we 
see from Fig. 3(a) that f\ > r c > r s , at least for 5 < 10. It 
means that a continuous transition to the chiral phase is 
not realized, and alternatively, a first-order transition to 
the helical phase occurs at a higher temperature. On 
the other hand, an inequality f\ < f c is satisfied for 
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FIG. 2: (Color online) The "transition temperatures" r c and 
f a as functions of v/u for the case of (a) u = 0.1, and (b) 
u = 0.5, where f c and f s are the para-to-chiral and the chiral- 
to-helical continuous transition points. 



v/u = 3.5. In this case, with decreasing the tempera- 
ture a continuous transition to the chiral phase occurs 
first, and then, the system goes into the helical phase 
through a first-order transition at f = f\ . 

The phase diagram of the 3D chiral GL model in the 
(f,v/u) plane is shown in Fig 0] for lager u. The dot- 
ted and the dashed curves represent continuous transi- 
tions which would occur if we would ignore the possi- 
bility of a first-order transition. The chiral phase pre- 
dicted by Onoda et al occupies the region between these 
two curves The solid curve represents the first-order 
transition from the paramagnetic phase to the helical 
phase which is newly found in this work. Although the 
range of the chiral phase is largely reduced due to the 
first-order transition, the chiral phase still persists for 
v/u > 2.5 in the case oiu= 200/(2tt 3 ) - 3.22. 

Thus, within the variational approximation, the chiral 
phase appears for sufficiently large u and v/u. The range 
of the chiral phase becomes wider for larger u, as can be 
seen from Fig[3J 

For the case of the XY spin (n = 2), the same analy- 
sis as was done in the Heisenberg case can also be per- 
formed by simply removing the r± term from the vari- 
ational Hamiltonian ^ and neglecting the 7 term: See 
appendix [XJ The results are qualitatively the same as 
those of the Heisenberg case, only the values of the tran- 
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FIG. 3: (Color online) The "transition temperatures" r c ,r s 
and fi as functions of u for the case of (a) v/u = 2.0, and 
(b) v/u — 3.5, where f c and r s are the para-to-chiral and 
the chiral-to-helical continuous transition points which would 
occur if the possibility of a first-order transition would be 
neglected, and f\ is the para-to-helical (or chiral-to-helical) 
first-order transition point. 

sition temperatures being different. Hence, we conclude 
that the chiral phase exists within the variational approx- 
imation for sufficiently large u and v/u even for the XY 
case. 



B. The u — > 00 limit 

In this subsection, we consider the behavior of the chi- 
ral GL model in the limit of u — > 00 with keeping r/u 
and v/u constant. In this limit, the model described by 
the Hamiltonian ([T]) reduces to the NLcr model [26| given 
by, 

?W = ^ J dr {{Va) 2 + (Vb) 2 + R[V(a X 6)] 2 }, 

(22) 

with conditions 

|a| 2 = |£| 2 = l, a±b. (23) 

From a simple saddle point analysis, one can derive the 
relation between the variables T and R of the NLcr model 
and the parameters of the chiral GL model. Indeed, the 
chiral GL model Hamiltonian ([!]) can be rewritten by 
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FIG. 4: (Color online) Phase diagram of the Heisenberg (n = 3) chiral GL model in the (f,v/u) plane for u = 200/(27r 3 ). 
The right figure is an enlarged view of the vicinity of the chiral phase. The red hatched (blue filled) area represents the chiral 
(helical) phase. The dotted and dashed curves represent continuous transition lines which would occur if the possibility of a 
first-order transition would be neglected in the analysis. The chiral phase predicted by Onoda et al occupies the region between 
these two curves fiol]. The solid curve represents the first-order transition line from the paramagnetic phase to the helical phase 
which is newly found in this work. The stability range of the chiral phase is largely reduced compared with that reported by 
Onoda et al. 



introducing r* — r/u and v* = v/u as 
U o = lj dr{(Va) 2 + (Vb) 2 } 

,2 ,.£2 



- / dr<, r | <r 



i 2 + P 



i * 

+ v 



(a-tj -a 2 b 2 



(24) 



Note that the second term of the r.h.s. becomes much 
larger compared with the first term in the limit of u — > oo . 
Therefore, one can apply the saddle point approximation 
to the partition function associated with (|2~4"1) . 



Va{r) / Vb(r)ex.p(-n c ). 



(25) 



Since the second term of (|24l) contains only on-site inter- 
actions, one can easily evaluate the minimization condi- 
tions for the second term as 



\b\ 2 



a Lb. 



(26) 
(27) 



In the limit of u — » oo, the functional integral of the 
partition function is approximated by an integral within 
a subspace constrained by the above conditions, yielding 



Z c ~ J J Va(r)Vb{r] 

|S| 2 =|?| 2 =1,5-L6 

xexpl— ^( d r 



2v* 



(Va) 2 + (Vbf 



(28) 



where we rescale the fields a(r) and b(r) so that they 
satisfy \a\ 2 — \b\ 2 — 1. Thus, one sees that the chiral 
GL model in the u — > oo limit reduces to the NLcr model 
with the correspondence 



T 8u-2w' 
R = 0. 



(29) 
(30) 



It also suggests that, in the limit of u — > oo, the proper- 
ties of the chiral GL model is independent of the value 
of v/u, depending only on the scaled parameter T given 
by Eq.fHU). 

Note that the chiral GL model corresponds to the 
R = sector of the NLcr model. Although David et 
al predicted that the stable chiral phase existed for suffi- 
ciently large R |2a|, there has been no report of the chiral 
phase for R = [27H30f . It indicates that for sufficiently 
large u there is very little chance for the chiral phase to 
be stabilized. 

By combining this result with that obtained from our 
variational analysis in the previous subsection, it seems 
most natural to expect that the chiral phase has the high- 
est chance to be stabilized for intermediate values of u 
and for larger values of v/u. Hence, in the following sec- 
tion, we perform extensive MC simulations focusing on 
such a parameter region in search for the possible chiral 
phase. 



IV. MONTE CARLO SIMULATION 
A. Method 

In this and following sections, we investigate the order- 
ing properties of the chiral GL model numerically by MC 
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simulations. For this purpose, the model is discretized on 
a 3D simple cubic lattice with lattice constant e as 



^4eWe(^ 



E 



bi 



+ r(a 2 + ty+u(a 2 + ty +v (a; • 6;) - 



(31) 



where Oj and bi are n-component vectors at the site i, 
while i + fi represents the nearest-neighbor site of i in 
the /i-direction (ju — x,y,z). We introduce rescaled pa- 
rameters as 



r' = e 2 r 


(32) 


u = eu 


(33) 


v = ev 


(34) 


_u i 
fl; = e 2 a, 


(35) 


&• - 


(36) 



The e dependence of the Hamiltonian (|3"Tj) is then re- 
moved as 

+ r> [(35) 2 + ^) 2 ]+«' [(«:) 2 + (^) 2 ] 2 

(^^-(^Wli (37) 



Such discretization and scaling procedure just corre- 
sponds to the cutoff procedure made for the continuum 
model (1). Setting A = ir/e, the parameters (r',u',v') in 
eqs. (|3"2"|) - (lMl) can be related to (f, u, v) in eg. (ITU as 



r O 



u O 



2^' 



v 



Note that a proportionality coefficient f/(27r 3 ) ~ 0.016 
means that u = 1 corresponds to a rather large value of 
u' ~ 62. 

We perform extensive MC simulations on the lattice 
chiral GL model described by (|3"7|) for both cases of the 
Heisenberg spin (n — 3) and the XY spin (n = 2) by 
using the standard Metropolis method. We consider r' 
as the "temperature" and the simulation is performed 
with the statistical weight of exp(— %) at each "temper- 
ature". The lattice is a L x L x L simple cubic lattice 
with 8 < L < 60. Periodic boundary conditions are im- 
posed in all directions. In updating a! i and b\ vectors, 
we adopt the polar coordinate in spin space and apply 
the type of Metropolis updating where appropriate win- 
dows are set for the proposed values of the variables so 
that the acceptance ratio becomes 0.25 ~ 0.6. When 



we simulate extreme cases where v! , v' and r' are much 
larger than unity, these adjustment procedures based on 
appropriate polar coordinate turn out to be efficient for 
thermalization. In case of the XY spin (n = 2), we also 
try the exchange of the directions of the and b\ vectors 
according to the Metropolis rule, which turns out to be 
efficient in relaxing the chirality vector aj x bi. 

Typically, our single MC run contains 6 x 10 6 MC steps 
per spin (MCS) at each temperature r' . In most cases, a 
gradual cooling protocol is employed. For large systems 
and near the transition temperature, we perform longer 
runs of up to ~ 1.8 x 10 7 MCS. In calculating physical 
quantities, initial (1 ~ 3) x 10 6 MCS are discarded for 
thermalization, and averages are calculated over subse- 
quent (3 ~ 5) x 10 6 MCS. Error bars are estimated by 
making 3^5 independent runs at each temperature and 
lattice size. 



B. Physical quantities 

In this subsection, we introduce various physical quan- 
tities measured in our MC simulation. 

We define the "specific heat" as the variance of the 
energy per spin, 



C^^({U 2 )-{U) 2 ) 



(39) 



Note that in the present model the specific heat de- 
fined by ([39]) is not equivalent to the one defined by 
the "temperature" derivative of the energy, 44, where 
E = jr < Ti > is the internal energy per spin, although 
both quantities are expected to exhibit similar singular 
behaviors at the transition. 

In order to measure the spin order, we define the spin 
order parameter by 



(38) We also define the chiral order parameter by 



^ W (( E«i x6 1) 2 > 



(40) 



(41) 



The finite-size correlation lengths of the spin and of the 
chirality are defined on the basis of the Ornstein-Zcrnikc 
form of the correlation function by 



1 



/ C a (0) 



- 1, 



(42) 



where a stands for either the spin (s) or the chirality (c), 
C a being the Fourier transform of the spatial correlation 
function, 



C s (q)^(K\ 2 + \b' q \ 2 ), 
C c (q)^(\(d'xb%\ 2 ), 



(43) 
(44) 
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where a' q , b' q and (a' x b') q represent the Fourier trans- 
form of a[, b[, and a - x & respectively. In our simulation, 
we take q m = (27r/L,0,0) corresponding to one of the 
minimum wavevectors compatible with periodic bound- 
ary conditions. 

We also measure the Binder ratios of the spin and of 
the chirality. For the spin, we have 



9s 



D 



D 



((E^) 4 > , ((Ez^) 4 ) 



k <(Ei«i) 2 ) 2 ((E» K) 2 Y 

while, for the chirality, we have 



9c 



D 



<(E 4 



(45) 



(46) 



The coefficients D[ a ^ and -D^ (a = s, c) are determined 
so that, in the thermodynamic limit, g a vanishes in the 
high-temperature phase and gives unity in the ordered 
phase. In the Hcisenberg case (n — 3), one has D — 
5/2, D ( 2 s) = 3/4, D[ c} = 5/2 and D { 2 c) = 3/2, while 



for the XY case (n = 2), one has d[ s) = 2, D y 2 > = 1 



Z?^ c) =3/2 and D 2 C > = 1/2 



(c) _ 



V. SIMULATION RESULTS: HEISENBERG 
(n = 3) CASE 

In this section, we present the results of our MC sim- 
ulation for the Heisenberg case (n = 3). As a typical ex- 
ample, we deal with the case of v! = 200 and v' /u' = 3.5 
here. For these parameters, the variational calculation of 
section IIIII predicts that the chiral phase is stabilized in 
a relatively wide temperature range between r' s ~ —23.6 
and r' c ~ —15.5. The relative difference between r' c and 
r s , a measure of the width of the chiral phase, might be 
given by 



Sr' 



\r' c + r' s \ 



(47) 



The variational calculation predicts a rather large value 
of Sr' ~ 0.42 for these parameter values. 

First, we show in Figj5] the temperature (r') depen- 
dence of the spin and the chiral order parameters. On 
decreasing the temperature, both order parameters rise 
up sharply around r' ~ —69.6, implying the occurrence 
of a phase transition around this temperature as can be 
seen from the figure. Each order parameter rises up at 
mutually close temperatures such that their relative dif- 
ference is much smaller than the one predicted from the 
variational calculation. 

Fig [6] exhibits the temperature (r') dependence of the 
specific heat. It exhibits only a single peak with no evi- 
dence of successive transitions. Note that the peak height 
of our largest size L = 60 is much larger than that of 




u' = 200 
v'/u' = 3.5 



FIG. 5: (Color online) The temperature (r') dependence of 
the spin and the chiral order parameters of the Heisenberg 
(n = 3) chiral GL model for various sizes. The Hamiltonian 
parameters are set v! = 200 and v'/u' = 3.5. 



140 
120 
100 
80 
60 
40 
20 




Heisenberg 

L=10 - 
1=16 -a- 
1=20 MM 
1=30 MM 
1=40 mm 
L=60 mm 




FIG. 6: (Color online) The temperature (r') dependence of 
the specific heat of the Heisenberg (n = 3) chiral GL model. 
The Hamiltonian parameters are set u = 200 and v'/u' = 3.5. 



other sizes. This may be a signature of a weak first- 
order transition occurring in the thermodynamic limit. 
In fact, as shown in figlZl the energy distributions for 
L = 60 shows double peaks characteristic of a first-order 
transition at the peak temperature, although the energy 
distributions of smaller sizes show only a single peak. 
In view of the fact that the parameter value studied 
here v'/u' = 3.5 is close to the mean-field tricritical line 
v'/u' — 4, the occurrence of a weak first-order transition 
seems consistent with the previous RG observation, since 
fluctuations extend the region of first-order transition. 

From the data of the order parameters and of the spe- 
cific heat, we now expect that, even if the chiral phase 



9 



Heisenberg 

L=20, r'=-69.50 

L=30, r'=-69.54 

L=40, r'=-69.56 

L=60, r'=-69.58 — 




19.7 



FIG. 7: (Color online) The energy distribution at the specific- 
heat peak temperature of the Heisenberg (n = 3) chiral GL 
model for various sizes. The Hamiltonian parameters are set 
u' = 200 and v' /u' = 3.5. 



exists, it is limited to a very narrow temperature range, 
never spreading as wide as Sr ~ 42% predicted from the 
variational theory. 

In order to further examine the possibility of an inter- 
mediate chiral phase, we show in FigsJS] and [S] the tem- 
perature (r') dependence of the Binder ratios and of the 
correlation-length ratios, respectively. As can be seen 
from FigJSl the Binder ratio of different sizes intersect 
almost at a common temperature for both cases of the 
spin and of the chirality at least for L < 40, whereas, for 
L = 60, the crossing point shows a downshift to lower 
temperature. This sudden change of the crossing behav- 
ior observed at L = 60 probably reflects the weak first- 
order nature of the transition observed in the specific heat 
and the energy distribution of this size. Indeed, both 
the spin and the chiral Binder ratios of L = 60 show a 
deep negative dip, which is a characteristic of a first-order 
transition. Recall here that just above a first-order tran- 
sition T = T+ the Binder ratio is expected to exhibit a 
divergent negative dip in the thermodynamic limit. The 
crossing temperatures of the correlation-length ratios de- 
pend on the system sizes only weakly for both cases of 
the spin and the chirality, as can be seen from FiglHl 

Although the difference in the crossing temperatures 
between different physical quantities, either the Binder 
ratio or the correlation-length ratio, or those between dif- 
ferent sizes, is sometimes of appreciable amount, we note 
that the difference in the crossing temperatures between 
the spin and the chirality is quite small for a given quan- 
tity and given sizes. In Fig. 10, we show the size depen- 
dence of the spin and the chiral crossing temperatures for 
both cases of the Binder ratio and the correlation-length 
ratio. The size is taken here as an average of the two sizes 
yielding the crossing point, L av — (Li + L 2 )/2. As can be 
seen from the inset, the difference is already very small 
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FIG. 8: (Color online) The temperature (r 1 ) dependence of 
the Binder ratio for the spin (a), and for the chirality (b), of 
the Heisenberg (n = 3) chiral GL model. The Hamiltonian 
parameters are set u 1 = 200 and v' /u' — 3.5. Insets are 
enlarged views of the transition region. 



even for small systems, of order 0.01 which corresponds 
to 0.1% relative difference, and tends to further decrease 
with increasing the system size. If we take the cross- 
ing temperatures of the Binder ratio between our two 
largest sizes L — 40 and L — 60, the spin and the chiral 
crossing temperatures, taken here as a measure of the 
respective transition temperature, are r' s = —69.599(2) 
and r' c = —69.598(2), which coincide within the errors. 
The relative difference between the spin and the chiral 
transition temperature is then limited to Sr' < 0.008%. 
Note that this upper limit is significantly smaller than 
the corresponding estimate obtained from the variational 
calculation Sr' ~ 42%. Hence, it turns out that the MC 
results are rather pessimistic about the occurrence of the 
chiral phase. 

We also try to estimate the bulk transition tempera- 
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FIG. 9: (Color online) The temperature (r') dependence of 
the correlation-length ratio for the spin (a), and for the chi- 
rality (b), of the Heisenberg (n = 3) chiral GL model. The 
Hamiltonian parameters are set v! = 200 and v' /u' = 3.5. 
Insets are enlarged views of the transition region. 



ture themselves by extrapolating the crossing tempera- 
tures to L — oo: See the main panel of Fig. 10. Such 
an extrapolation procedure is hampered somewhat by 
the sudden change of the behavior observed at L = 60. 
This effect seems relatively minor for the correlation- 
length ratio. By performing a power-law fit of the form 



XL) 



-XL 



oo) + cL to the data of £ S /L 
and of &./L, we get r' cross (L = oo) = -69.603 ± 0.006 
from £ S /L, and r' cross (L = oo) = -69.600 ± 0.004 from 
£ C /L. These estimate are consistent with our estimate 
based on the Binder ratio given above. 

We also perform similar MC simulations for other pa- 
rameter values, including v! = 1,10,100,200,1000 and 
v' /u' — 2.0, 3.5. For all these parameter values, there is 
no evidence of the chiral phase as in the case of vl = 200 
and v' /u' = 3.5 shown above. In case of it' = 1, 10 and 
v 1 ' jvl = 3.5, signatures of a strong first-order transition 
are observed. By contrast, for other parameter values, 
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FIG. 10: (Color online) The size dependence of the crossing 
temperatures of the Binder ratio and the correlation-length 
ratio for both cases of the spin and the chirality of the Heisen- 
berg (n — 3) chiral GL model. The size is taken here as a 
mean of the two sizes L av — (Li + La)/2, where L\ and L2 
are the two liner system sizes yielding a crossing point. The 
crossing points associated with the L = 60 data are indicated 
by arrows. Inset exhibits the size dependence of the differ- 
ence between the crossing temperatures of the spin and of the 
chirality for both cases of the Binder ratio and the correlation- 
length ratio. 
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TABLE I: (Color online) The spin and the chiral transition 
temperatures and its relative difference for various parameter 
values of the Heisenberg (n = 3) lattice GL Hamiltonian, 
estimated from the crossing points of the Binder ratio of the 
two sizes L = L\ and Li. 



the transition seems to be continuous, at least in the 
range of lattice sizes studied here. We summarize our 
estimates of the spin and the chiral transition temper- 
atures in table HI where each transition temperature is 
estimated from the crossing temperatures of the Binder 
ratio between the sizes L = L\ and L = L^- For all cases 
studied, the difference Sr' is less than 0.1%. 

Finally, we compare the transition temperature of the 
lattice chiral GL model as estimated from our MC with 
that of the NLct model. In FigfTTl we plot the transition 
temperature estimated from our MC versus the parame- 
ter 8u' — 2v'. For u' > 100 , the transition temperature 
can be well fitted by the NLcr model relation (|29p , if one 
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FIG. 11: (Color online) The transition temperatures of the 
Heisenberg (n = 3) chiral GL model estimated from our 
Monte Carlo simulations are plotted versus the GL Hamil- 
tonian parameter 8u' — 2v' . The line shows the relation (|29p 
with T — 3.062, which is the transition temperature of the 
NL<r model reported in ref.[29l] (Note that we multiply the 
transition temperature by two due to the difference in the 
definition). 



identifies T in Eq. (j2U)l as the transition temperature of 
the lattice NLct model reported in ref.[29|]. It also indi- 
cates that the spin and the chirality order simultaneously 
and there exists only single transition in the large 8u'—2v' 
region, since a common belief is that there is only single 
transition in the NLer model with R = [26T4291 . 



VI. SIMULATION RESULTS: XY (n = 2) CASE 

Next, we present the result of our MC simulation for 
the XFcase (n = 2). As a typical example, we deal with 
the case of u' — 200 and v'/u' — 3.5 again. For these 
parameters, the variational calculation of section lHTl (and 
appendix [5} predicts that the chiral phase is stabilized 
in a wide temperature range between r' s ~ —21.6 and 
r' c ~ —13.4. The relative difference between r' c and r' s is 
Sr ~ 0.47. 

As we discussed in section IIII1 the model reduces to 
the NLcr model in the limit of u' — > oo. Previous stud- 
ies showed that, for the case of n — 2, the NLct model 
discretized on a 3D lattice exhibited a single first-order 
transition into the magnetic ordered state [53, ■ Mean- 
while, the variational calculation predicts the stable chi- 
ral phase for sufficiently large u' and v'/u'. 

First, we show the temperature (r') dependence of the 
spin and the chiral order parameters in Fig |12l On de- 
creasing the temperature, both order parameters rise up 
sharply around r ~ —47.7, implying the occurrence of 
a phase transition around this temperature. Each order 
parameter rises up at mutually close temperatures such 
that their relative difference is much smaller than the one 
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FIG. 12: (Color online) The temperature (r 1 ) dependence of 
the spin and the chiral order parameters of the XY (n — 2) 
chiral GL model for various sizes. The Hamiltonian parame- 
ters are set u' = 200 and v'/u' = 3.5. 



predicted from the variational calculation. The observed 
onset of the order parameters seems steeper than the one 
observed in the Heisenberg case. Such a sharp rise of the 
order parameters is suggestive of a first-order transition. 

We show in Fig[T2]the temperature (r') dependence of 
the specific heat, in Fig|l4]the energy distribution near 
the transition temperature. The specific heat exhibits 
only a single peak with no evidence of successive transi- 
tions. Note that the peak height grows rapidly with in- 
creasing the system size L, consistently with a first-order 
nature of the transition. In fact, the energy distribution 
near the transition point shows double peaks character- 
istic of a first-order transition. Due to the difficulty in 
thermalizing the system exhibiting a rather strong first- 
order transition, the lattice sizes in the XY case are re- 
stricted to be smaller than those in the Heisenberg case. 

Although we observe a rather strong first-order tran- 
sition from the paramagnetic phase into the helical mag- 
netic phase, there still remains a possibility that a chiral 
phase is stabilized between the tiny temperature region 
between the paramagnetic phase and the helical phase. 
Since the variational calculation of section UTTI predicted a 
continuous transition from the paramagnetic phase to the 
chiral phase, we examine here the possibility of a continu- 
ous transition occurring at a temperature higher than the 
first-order transition temperature. If such a transition re- 
ally occurs, chirality-related dimensionless quantities of 
various sizes, e.g., the chiral correlation-length ratio and 
the chiral Binder ratio, are expected to exhibit a cross- 
ing behavior at a higher temperature than the first-order 
transition temperature. 

Fig. [TBI exhibits the temperature (r') dependence of 
the spin and chiral correlation length ratios. One can see 
from the figure that, with decreasing the temperature, 
both £ S /L and £, C /L rise up sharply around r' ~ —47.75, 
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FIG. 13: (Color online) The temperature (r ) dependence of 
the specific heat of the XY (n — 2) chiral GL model. The 
Hamiltonian parameters are set u' = 200 and v' /v! = 3.5. 



5 
4.5 
4 
3.5 
g 3 



2 

1.5 
1 

0.5 






L 


XY 






L=8,r'=-47.55 






L=12,r'=-47.69 






L=16,r'=-47.74 






L=20,r'=-47.76 






. u = 200 . 






f \ v' Iv! = 3.5 






— 


-8.6 -8.4 


-8.2 -8 - 

E 


7.8 -7.6 -7.4 



FIG. 14: (Color online) The energy distribution near the 
specific-heat peak temperature of the XY (n — 2) chiral GL 
model for various sizes. The Hamiltonian parameters are set 
u = 200 and v'/u' = 3.5. 



close to the specific-heat peak temperature. The cross- 
ing temperatures between our two largest sizes L = 16 
and L = 20 are r' s = —47.758(5) for the spin and 
r' c — —47.753(3) for the chirality. The relative difference 
between the spin and the chiral crossing temperature is 
Sr' < 0.03%. The spin and the chiral crossing points co- 
incide within the errors, and are close to the first-order 
transition temperature estimated above. 

We show in Fig. Q1)]the temperature (r') dependence 
of the spin and the chiral Binder ratios. As expected for 
a first-order transition, the Binder ratios exhibit a deep 
negative dip, which grows with L. In addition, the Binder 
ratios exhibit a crossing on the positive side of g at a 
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FIG. 15: (Color online) The temperature (r') dependence of 
the correlation length ratios for the spin (a), and for the chi- 
rality (b), of the XY (n = 2) chiral GL model. The Hamilto- 
nian parameters are set u — 200 and v'/u = 3.5. Insets are 
enlarged views of the transition region. 



temperature slightly below the dip temperature: See the 
insets. The observed behavior of the chiral Binder ratio 
is hardly compatible with a continuous chiral transition 
occurring at a temperature higher than the first-order 
transition temperature in the thermodynamic limit. The 
crossing temperatures between our two largest sizes L = 
16 and L = 20 are r' s = -47.777(9) for the spin and 
r' c = —47.776(8) for the chirality. The relative difference 
between the spin and the chiral transition temperatures 
is then estimated to be Sr' < 0.04% 

We also perform similar MC simulations for the case 
of v'/u' = 2.0 and u' = 100, located rather far from the 
mean-field tricritical line v'/u' — 4.0. Even in this case, 
we find a signature of a first-order transition such as the 
double-peak structure in the energy distribution. Fur- 
thermore, no evidence of successive transitions are ob- 
served. The behaviors of the correlation-length ratios 
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FIG. 16: (Color online) The temperature (r') dependence of 
the Binder ratio for the spin (a), and for the chirality (b), of 
the XY (n — 2) chiral GL model. The Hamiltonian param- 
eters are set it' = 200 and v'/u' = 3.5. Insets are enlarged 
views of the transition region. 

and the Binder ratios are qualitatively the same as in 
the case of v'/u' = 3.5 and u' — 200, and the relative 
difference between the spin and the chiral transition tem- 
perature is limited to Sr < 0.05%. 

VII. SUMMARY AND DISCUSSION 

In this paper, the nature of the phase transition of reg- 
ularly frustrated vector spin systems in three dimension 
was investigated based on the chiral GL model both by 
analytical calculations and Monte Carlo simulations. 

We first performed a variational calculation which was 
an extension of the previous calculation by Onoda et 
a?[Tfj|. and showed that the chiral phase was stabilized in 
a certain restricted parameter range. We considered the 
possibility of a first-order transition, which was not con- 



sidered by Onoda et al. We then found that a first-order 
transition indeed occurred in this model significantly re- 
ducing the stability range of the chiral phase. Yet, we 
observed that the chiral phase still persisted for suffi- 
ciently large u and v/u within the variational approxi- 
mation. We also showed that in the limit of u —> oo 
the chiral GL model reduced to the NLcr model without 
the coupling between the chiralities (R = in eq.(22)). 
Previous analysis of the NLcr model predicted only a sin- 
gle phase transition directly from the para to the helical 
phase, which means that there is very little chance for 
the chiral phase to be stabilized for sufficiently large u. 

With reference to these analytic results, we next per- 
formed extensive MC simulations on the lattice dis- 
cretized version of the 3D chiral GL model in search for 
the possible chiral phase. In contrast to the expecta- 
tion based on the variational results, however, we did 
not find any numerical evidence of the chiral phase for 
both cases of the Heisenberg model (n = 3) and the XY 
model (n = 2). From the data of the correlation length 
ratios and the Binder ratios, we conclude that for all cases 
studied the stability range of the chiral phase, if any, is 
less than 0.1% in the relative temperature width. 

A possible cause of the appearance of the chiral phase 
in the variational calculation might be the fact that, for 
the chirality, only the linear term 5a x Sb y — 5a y db x is 
contained in the variational Hamiltonian ([5]), while the 
quadratic terms are also present for the spin. This imbal- 
ance inherent to the variational calculation might lead to 
an underestimate of chirality fluctuations compared with 
spin fluctuations. 

Although we could not find any evidence of the chiral 
phase in the chiral GL model (JTJ, there still remains a 
possibility of the chiral phase originated from some other 
mechanisms not included in the chiral GL model. If we 
consider, for example, the direct interaction between the 
chiralities, i?[V(a x b)] 2 , it enhances the ordering of the 
chirality. For sufficiently large R, the chiral phase is triv- 
ially stabilized. In fact, David and Jolicoeur showed on 
the basis of the NL a model that the chiral phase is sta- 
bilized for sufficiently large R [26[. However, for smaller 
R, which is more realistic, their analysis indicated that 
the chiral phase disappeared and there was a single tran- 
sition from paramagnetic phase to the helical phase ptjj . 
To get further insight into the effect of the i?-term, we 
also performed a MC simulation of the chiral GL model 
with a weak chirality-chirality interaction (the i?-term) 
with u' = 200, v'/u' = 2.0 and R = 0.1. The results turn 
out to be qualitatively the same as those of the the origi- 
nal chiral GL model with R = 0, and Sr' estimated from 
the crossing temperature of the Binder ratios is also very 
small, Sr' < 0.02%. 

Another mechanism to stabilize the chiral phase in 
3D regularly frustrated system was proposed by Villain 
[3ll ]. He suggested that a chiral phase might be real- 
ized in quasi-one-dimensional XY spin systems when the 
inter-chain coupling is sufficiently weak |3l| . In purely 
one-dimensional frustrated XY spin systems, it has been 
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known that, with decreasing the temperature, the chiral- 
ity correlation length increases exponentially, while the 
spin correlation length diverges as a power law. By taking 
into account the effect of the weak inter-chain coupling, 
Villain predicted that the chiral phase could exist in an 
intermediate temperature regime above the helical phase. 
As far as the authors know, however, a direct numerical 
evidence of such a chiral phase in quasi-one-dimensional 
system is still lacking. Thus, it is an interesting open 
problem to observe the chiral phase by numerical simula- 
tions of such quasi-one-dimensional frustrated spin mod- 
els. 
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Appendix A: Details of a variational approximation 
1. Free energy 

In this Appendix, we explicitly show the form of the 
trial free energy Fq + {% — V.o)o- 

Since the variational Hamiltonian (6) is diagonalized 
with respected to a, /3 and 7 as given in (8), the term J-"o 
is easily calculated as 

To = - log I J Y[ d5a q dSb q exp(-'Ho) j 



V(<Z 2 + r,| - h K /2) 



+ log 



V(q 2 + r,| + h K /2) 



log 



W + rj.) 

(Al) 



In calculating the second term (H — "Ho)o, we rewrite 
eq.(l) in terms of a, /3 and 7 defined by eq.(7), to get 
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H = I f dxy (a 2 + /3 2 + f ) + (Va) 2 + (V/3) 2 + (V7) 2 + (u - ~«) (a 4 + /3 4 ) + 2 ^ + i« ) a 2 /3 2 

+W 4 +2ln- i« J f (a 2 + /3 2 ) + v\2 lll2 {d ■ 0) + ( 7l 2 - ji){a 2 /3i - ai#j)] L. (A2) 



r 



By substituting a = Sa + A, etc, into (A. 2), (T-L — Ho)o is calculated via simple Gaussian integrals as 
I 



r - ( <1 ~ y ) + 4 (« - I) A 2 + 2 (u + \ ) if + 2 { „ - - ) C* 



(-3 



A 2 +4u--B 2 + 2 



(-1) 



r — r±_ 



2 ( u -j) {A 2 + BA + AuC 2 ] • 57-a)o 



-2[u-- 



+r(A 2 + B 2 + C 2 ) + (u - \ (A 4 + B A ) + 2 (u + J A 2 ^ 2 + tiC 4 



-2 (« - |« ) C 2 (A 2 + £ 2 ) +w 



2C 1 C 2 (A • B) + (C? - Cf )(A 2j B x - A X S 2 ) 



r 



(A3) 



In the case of Ai = m, A 2 = 0, 5 = 0, C = 0, as assumed in §IIIA, the trial free energy is simplified as 
I 



VA 3 f 1 { 

T Q + (H- n Q )o = --^2- J Q M lo e 

VA 3 



1 



q 2 +f|| - fr K /2 



log 



1 



<7 2 + f|| + ft„/2 



log 



1 



4tt 2 



r — f|i H — — + 4 I £t — 7 ) m 



-2 



f — f 11 — + 2 u + — m 



s,2 



r — f j_ + 2 I u I m 



+2w(7 2 + 2 ( u - ^ ) cr 2 (cr 2 + ct^) + fra 2 + ( u - -€) to 4 I + const., (A4) 



where we introduced the cutoff wavevector A, and scaled various parameters as in pip . By taking the 
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derivatives of F$ + (H — Ho)o with respect to f | , f±,h K 
and to, and setting them to zero, we can get the condi- 
tions for the optimal parameter values given in ([T2|) - ([T5l) . 



2. Relation between f a and f c 

In this Appendix, the relation between the para-to- 
chiral and the chiral-to-helical continuous transition tem- 
peratures, f c and f s , is investigated. We prove here the 
existence of a critical value u c such that f s > f c for u < u c 
and f s < f c for u > u c . 

To simplify the notation, we define a function 



9( r ) = I ~T~ d 1 
Jo T + r 



1 — -y^a-tan 



1 



(A5) 



By using <?(r), the variances of a, /3 and 7 are calcu- 
lated as 



= 2 5 (f|| - fe«/2), 
aj=2g(f ]] +K/2), 



ol = 2g(f ± ). 



(A6) 



The conditions of optimal parameter values eqs. (fT12)) - 
14]) may be given with to = by 



n : ' - '- } ( : '>« - 1 ) ff(r|| - K/2) + g(f {l + h K /2) 

(A7) 

g(f l{ - K/2) + g(f || + h K /2) 
+ 8%(f x ), (A8) 
1 = (35 - 4fi) ( ff (f|| - h K /2) - g(f {l + K/2)) IK- 



f± = r + 4 ( u 



(A9) 



The equations determining the critical parameters 



~( c ) 



and f,| , eqs.(20) and (21), are also given by 



(3«-4u)/(f[| C) ) = 1 



(3i;-4fi) l 3 (0)- 5 (2r~[ | s) ) 



2f 



(A10) 
(All) 



where the function f(r) is defined by 



er„\- d 9 
f(r) = --T 
dr 



111 
:atan- 



y/r y/r 1 



(A12) 



Originally, eqs. (|A7p - (|A9I) are the equations determin- 
ing the variational parameters f || , fj_ and /i K as functions 
of the Hamiltonian parameters f, u and v. One can also 
regard f given as a function of h K , u and v. In order 
to discuss the relation between f s and f c , we examine 
here the behavior of f(h K ; u, v) around h K = for given 
u and -5. Because of symmetry, df/dh K is equal to zero at 
h K = (at f = f c ). Hence, the sign of the second deriva- 
tive of f(/i K ) at the chiral transition point h K — >• deter- 
mines the relation between f s and f c . Namely, f s < f c if 
d 2 f/dh 2 K < 0, and f s > f c if d 2 f/dh 2 K > 0. 

By taking the second derivatives of equations (|A7[) - 
(|A9|) with respect to h K and setting f|| = f^\ f± = f^ 
and h K = 0, the second derivative d 2 f/dK K is given by 



d 2 f 



l + (16u-2u)/(ff; c) ) 



r(fj c) ) 

12/'(f[, c) ) 



+ 1 + (4S + C) /(f) /(r " } ' (A13) 

where f'(r) and /"(r) are the first and the second deriva- 
tives of f(r), respectively. By use of eq. (|A10l) . eq. (IA13|) 
can be rewritten as 



Cr(s) 



/'(f) 



dhl C 2 (s) ftfftf,,^)) 



C 2 (s) 



r(rj e) ) 

>(rj c) ) 



(A14) 

where we set s = u/u, while the functions Ci(s) and 
C*2(s) are defined by 



c 2 ( s ) 



16- 2s 



3s 
4- 



3s -4 



4s 



(A15) 



(A16) 



Note that Ci (s) and C2 (s) are positive for 4/3 < s < 4. 
It is easily confirmed that Ci(s)/C2{s) is an increasing 
function of s, lying in the range 1/3 < C\(s)/C<2{s) < 2/3 
for 4/3 < s < 4. One can also show that the function 
[/'( r )] 2 /[/( r )/"( r )] i s an increasing function of r in the 
range between 1/3 and 2/3 for r > 0. Hence, there exists 
a critical value r*(s) such that d 2 f/dh 2 K > for fy < 

r*(s), and d 2 f/dh 2 < for fjj c ' > r*(s), where r*(s) is 
given by the solution of the equation 



Cx{s) [f(r*)Y 



0. 



(A17) 



C 2 (s) f(r*)f"(r*) 

We have r c > f s if and only if fjj > r*(s), while 
f c < f s if < r*(s). 



17 



From eq. (|A10l) . fy is an increasing function of u for 

fixed s = v/u, the critical value u c (s) corresponding to 
r* is given by 



section Hill the variational Hamiltonian for the XT case 
may be given by 



(3s-4)/(r*)- 



(A18) 



Thus, the same conclusion can be restated in terms of 
u as r c > f s if and only if u > u c (s). This critical value 
u c (s) is an increasing function of s = v/u and is in the 
range < u c < oo for 4/3 < v/u < 4. 



Hq 



V 



n 



r —J Sd q ■ Sa-q 



+ [q' + r\\+-£)6l3 q -6P- q . (A19) 



3. Variational approximation for the JfKspin 
(n = 2) 

In this Appendix, we show the results of the variational 
calculation for the XT case (n = 2). As mentioned in 



As in the Heisenberg case, we assume that A\ = 
m > 0, A2 = and B = 0. The trial free energy 
Jo + (H — Ho)o is then calculated as 



VA 3 
2tt 2 



dql log 



1 



q 2 + f|, - h K /2 



log 



q 2 + f|| + h„/2 



yA 3 
'4tt 2 ' 

+2 ( u 



•i' 



h 4 {t m 



4 

2 (u 



f — f|| ■ + 2 ( U + — ) TO 



2 2 , 2 

o a <Jp + rm 
I 



-v I to 4 > + const. 



(A20) 



By taking the derivatives of Tq + (H — Ho)o with re- 
spected to r 11 , h K and to, and setting them to zero, we 
get the following conditions for the optimal parameters 
values, 



3u-~) (m 2 + a 2 a + af 




(A21) 
(A22) 

(A23) 



The difference between these equations and the cor- 
responding equations of the Heisenberg case (fT2"]) - (fT5)) is 
only the terms related to r± or a 2 which are simply ab- 
sent in the XY case. Especially, eq. (|A22p is exactly the 
same as eq. (fl"4")) . Therefore, when v/u > 4/3, there is a 
continuous transition to the chiral phase at the same crit- 
ical value of f || = f|| as in the Heisenberg case. The con- 
tinuous transition between the chiral phase and the heli- 

( s) 

cal phase also occurs at the same critical value fy = r\ 
as in the Heisenberg case. On the other hand, if one an- 



alyzes how the ordering proceeds when the temperature 
f is varied, the transition temperatures of the XY spins 
differ from those of the Heisenberg model quantitatively. 
Both r s and r c are generally higher than the correspond- 
ing transition temperatures of the Heisenberg model. 

By use of the same technique as in Appendix IA 21 one 
can show that there is the critical value u c such that 
f c > f s for u > u c . From eqs. (IA21l) and (|A22|) . the same 
form of equation for df 2 /dh 2 K with (|A14[) can be derived, 
where the functions C\{s) and C2(s) are different from 
those of the Heisenberg case and are given by 



cm 



C 2 (s) 



1 

12 



12 



3s -4 



3s 



(A24) 
(A25) 



Again, Ci(s) and 62(5) are positive for 4/3 < s < 4, 
while the ratio C\{s) / C2{s) has the same characteristics 
as those of the Heisenberg case such that it is an increas- 
ing function of s satisfying 1/3 < Ci(s)/C 2 (s) < 2/3. 
Therefore, we reach the same conclusion as in the Heisen- 
berg case, i.e., r c > r s , if and only if u > u c (s). Note 
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that the value of u c (s) differs from that of the Heisenberg 
model due to the difference in the function C\{s)/C-i{s). 

In the same manner as in the Heisenberg case, we find 
that a first-order transition to the helical magnetic phase 
occurs and this first-order transition reduces the stabil- 
ity range of the chiral phase. However, even in case of 



the XY spin, the chiral phase still persists for a certain 
parameter range. In Fig |17l we show a typical phase di- 
agram of the XY case. One can see qualitatively similar 
structure to the one of the Heisenberg spins (see also Fig. 
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FIG. 17: (Color online) Phase diagram of the XY (n = 2) chiral GL model in the (r,v/u) plane for u = 200/(2tt 3 ). The right 
figure is an enlarged view of the vicinity of the chiral phase. The red hatched (blue filled) area represents the chiral (helical) 
phase. The dotted and dashed curves represent continuous transition lines which would occur if the possibility of a first-order 
transition would be neglected in the analysis. The chiral phase predicted by Onoda el al occupies the region between these two 
curves[l(J. The solid curves represents the first-order transition line from the paramagnetic phase (or the chiral phase) to the 
helical phase. As in the Heisenberg case, the stability range of the chiral phase is largely reduced compared with that reported 
by Onoda et al. 



